Anomalous Scaling and Fractal Dimensions 
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The relation between critical exponents, characterizing a continuous phase transition, and the 
fractal structure of physical lines, proliferating at the critical point, is established by considering the 
two-dimensional O(N) spin model for which many exact results are available. 



A field theory describes a grand canonical ensemble 
of fluctuation loops of arbitrary shape and length, a so- 
called loop gas. In the context of quantum field theories, 
the fluctuating loops physically represent the worldlincs 
traced out by the quantum particles, while in a classical 
context they represent, for example, high-temperature 
(HT) graphs (in the case of spin models), current loops 
(in the case of the Ginzburg-Landau theory) or closed 
stream lines (in the case of the complex \(f>\ 4 theory). In 
a recent comment, Prokof'ev and Svistunov Q pointed 
out that a relation proposed by Hove, Mo, and Sudb0 
between the anomalous scaling dimension of a field 
and the fractal or Hausdorff dimension of critical loops is 
incorrect. They concluded furthermore that the anoma- 
lous scaling dimension cannot be deduced from simula- 
tions of closed loops only. They supported their criticism 
by Monte Carlo simulations of the HT representation of 
the three-dimensional complex |(/>| 4 theory and reported 
the value D = 1.7655(20) as fractal dimension of the HT 
graphs at the critical point. 

As was first pointed out by Helfrich and Miiller 3] , the 
HT expansion of an AT- vector spin or O(N) lattice model, 
which can be visualized by closed graphs along the bonds 
on the underlying lattice 0, describes at the same time 
a loop gas of sterically interacting physical lines. Specif- 
ically, the HT expansion of the three-dimensional XY 
model (N — 2), which belongs to the same universal- 
ity class as the \4>\ 4 theory, simultaneously represents the 
vortex loop gas of a bulk superconductor [5j . As the mag- 
netic interaction is screened, vortices in superconductors 
have short-range interactions that can be accurately de- 
scribed by a steric repulsion. The equivalence implies 
in particular that the fractal dimensions of HT graphs of 
the three-dimensional XY model and vortex loops in bulk 
superconductors, studied in Ref. 0], coincide. Prokof'ev 
and Svistunov |lj therefore rightly compare their results 
with those of Ref. 0, where the value D = 1.92(4) was 
reported for magnetic vortex loops at the critical point. 

In their reply to the comment pj, Hove and Sudb0 Q 
questioned the validity of the criticism leveled at their 
earlier work. By considering the HT representation of 
the two-dimensional O(N) spin model for — 2 < N < 2, 
for which many exact results are available, we give in 
this short note analytic arguments supporting the find- 
ings of Prokof'ev and Svistunov and, hopefully, settle the 



issue. The HT graphs of the two-dimensional Ising model 
(N = 1) physically represent Peierls domain walls, sep- 
arating spin clusters of opposite orientation on the dual 
lattice. These closed lines, which we recently investigated 
numerically in Ref. Q, form the two-dimensional coun- 
terpart of magnetic vortex loops in three dimensions 

A loop gas can be conveniently characterized by the 
loop distribution £ n , giving the average number of loops 
of n steps per unit volume. Close to the critical point 
K c , £ n asymptotically takes a form similar to the cluster 
distribution in percolation theory [8j, 
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Here, 6 is the line tension, r and a are two exponents 
characterizing the distribution, and K is the tuning pa- 
rameter. When the line tension is finite, the Boltzmann 
factor in the distribution Q exponentially suppresses 
large loops. Upon approaching the critical point, 9 van- 
ishes at a rate determined by the exponent a. At K c , the 
loops proliferate as they can now become arbitrary long 
without energy penalty. The remaining factor in the loop 
distribution is an entropy factor, giving a measure of the 
number of ways a loop of n steps can be embedded in the 
lattice. The entropy exponent r is related to the fractal 
dimension D of the loops at the critical point via 



t = d/D + 1, 



(2) 



where d denotes the dimension of space (in the case of 
classical theories) or spacetime (in the case of quantum 
theories) . The entropy factor decreases with increasing n 
because it becomes increasingly more difficult for paths 
to close the longer they are. 

In the limit N -> 0, the HT graphs of the 0(A) 
model reduce to self-avoiding random walks (SAWs) 
which physically represent polymer rings in good solvents 
9]. Using this equivalence, de Gennes determined the 
fractal structure of SAWs in terms of the critical expo- 
nents of the 0(N — > 0) model. Generalizing his results 
to arbitrary — 2 < N < 2, we showed in a recent Letter 
[lf)j | that the fractal dimensions D of the HT graphs and 
the exponents a of the loop distribution JTJ are related 
to the correlation length exponent v of the O(N) model 
in the following way: 



= 1/o-D. 



(3) 
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TABLE I: Critical exponents of the two-dimensional critical 
O(N) spin models, with N = —2,0, 1,2, 00, respectively, to- 
gether with the the fractal dimension D of the HT graphs as 
well as the two exponents a and i9. 



FIG. 1: A SAW in two dimensions returning to a site adjacent 
to its starting point x. 



For SAWs (N — > 0) a = 1, but in general a takes a 
different value. 

In the context of SAWs it is well established that loops 
alone do not yield all the critical exponents of the uni- 
versality class defined by the 0(N — > 0) model [9j. To 
that end, also the total number z n = z„(x, x') of 
SAWs of n steps starting at x and ending at an arbi- 
trary site x' is needed. Because of translation symme- 
try z„(x, x') = 2 n (|x — x'|), and z n does not depend on 
x. The ratio of z„(x, x') and z n defines the probability 
P„(x, x') of finding a path connecting x and x' in n steps. 
As n — > 00, it scales as 01 



P n (x,x') = z n (x,x')/z„ ~ n- d ' D P f|x - x'l/n 17 



(4) 

with P a scaling function, while the number z n scales as 

(5) 



with 1? an exponent which is to be identified. Since the 
number of possible rooted open paths with no constraint 
on their endpoint increases with the number n of steps, 
$ is expected to be positive. 

The loop distribution Q is related to the number of 
paths £n(x, x + ai) = z n (a) of n steps returning to a site 
x + ai adjacent to the starting point x (see Fig. and 
Ref. 0) through 
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must vanish for t — > and behave for small argument t as 
P(t) « t with an exponent determined by the asymptotic 
behavior (JSJ of the number z n of open paths at the critical 
point (0 = 0). 

The fractal dimension D and the exponents a and $ 
determine the critical exponents of the theory. As for 
self-avoiding random walks, the relations can be derived 
by writing the correlation function G(x, x') in terms of 
z„(x,x') as G(x,x') ~ ^ n z„(x,x'), where because of 
translational invariance, G(x, x') = G(|x — x'|). When 
evaluated at the critical point, where G(x, x') ~ l/|x — 
x'| d ~ 2+I? , this gives the relation proposed by Prokof'ev 
and Svistunov [l| 
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with i any of the (positive or negative) directions on the 
lattice, V the lattice volume, and a the lattice spacing. 
The latter serves as a microscopic cutoff, whence z n (a) 
rather than z n (0) appears when closing open paths. The 
factor 1 jn in Eq. © accounts for the fact that a loop can 
be traced out starting at any lattice point along the loop. 
As first shown by McKenzie and Moore for self-avoiding 
random walks [13|, consistency of Eqs. © with z n (a) = 
z n P n {a) and Q requires that the scaling function P(t) with x 



Using high-precision data for 77, together with their esti- 
mate for D, they arrived at the estimate d — 0.1965(20) 
for the three-dimensional |0| 4 theory. The relation pro- 
posed by Hove, Mo, and Sudb0 corresponds to set- 
ting -& to zero in Eq. J2J, which in general is not al- 
lowed. Given the exact values for 77 Il4l| and the fractal 
dimensions D of the HT graphs |l5j, 1? can be deter- 
mined exactly for the two-dimensional O(N) model (see 
Table QJ. Through the exact enumeration and analysis 
of the number z n of self-avoiding walks on a square lat- 
tice up to length 71, the expected value d/D = 11/32 for 
N = has been established to high precision [16j . while 
our Monte Carlo simulation of the high-temperature 
representation of the two-dimensional Ising model un- 
ambiguously shows that t? is nonzero. In their numerical 
study of the three-dimensional |</>| 4 theory Prokof'ev 
and Svistunov considered not the asymptotic behavior 
of the number z n of open graphs to verify the value 
i? = 0.1965(20), but the behavior of the probability @ 
for small arguments. Both should give the same results, 
according to our analysis. 

Relation can, incidentally, be derived by using the 
second-moment definition of the correlation length £, 



£ x G(x) ' 



(8) 



|x|. Finally, using the definition of the sus- 
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ceptibility \, x — Sx' G(x,x'), which diverges as \ ~ 
|if — if c | -7 , we obtain 

7 = {D + $)/ctD. (9) 

This relation generalizes one due to Cloizeaux 0] f° r 
SAWs for which a = 1. The explicit expressions for ^,77, 
and 7 satisfy Fisher's scaling relation, 7/1/ = 2 — 77. The 
critical exponents are a function of the two independent 
variables z\ = D + $ and 22 = o~D which determine the 
anomalous scaling dimension of the <f> and (j) 2 fields: 

d^ = \{d-2-rj) = \{d- Zi), 

d^ = d =d — Z2, (10) 

respectively. 

With d set to zero, Eq. reduces to the relation 
7 = 1/cr proposed by Nguyen and Sudb0 in Ref. (lif . 
Given that $ is in general positive, that relation, too, 
is generally incorrect. It is, however, not impossible for 
i? to be zero. A class of models with the special value 
i? = is provided by a noninteracting Bose gas in d space 
dimensions with modified energy spectrum e(k) <x k D 
with D < d < 2D at fixed particle number density [l9| . 
This class belongs to the universality class of the spherical 
model, obtained by taking TV — > 00 in the O(N) model, 
with long-range interactions |20j |. The critical exponents 
of this universality class are exactly known [ilj . Those of 
the two-dimensional spherical model without long-range 
interactions, corresponding to D = 2 and zero critical 
temperature, are included in Table [I] 

The presence of the exponent t? is specific to having 
geometrical objects, viz. lines, that can be open or closed 
as it provides consistency when closing open lines. For 
geometrical objects where the notion of open and closed 
does not apply, such as for spin clusters, the analog of 
the exponent 1? is absent. This is, for example, the case 
in the Fortuin-Kasteleyn formulation of the Q-stete Potts 
model 22] , and for the improved estimators in the cluster 
update Monte Carlo algorithms |23| . 

In conclusion, using exact results for the two- 
dimensional O(N) model with — 2 < N < 2, we provided 
analytic support to the findings of Ref.~Q] that the full 
set of critical exponents for a given theory cannot be de- 
termined from the loop distribution alone. For the full 
set also the asymptotic behavior of the number of open 
paths is needed. Generalizing results known for SAWs, 
we showed that the exponent 1? governing this behavior is 
connected to that of the probability -P„(x, x') of finding a 



path connecting x and x' in n steps for x' — > x, i.e., in the 
limit of closing open paths. Finally, we showed that this 
connection, together with a relation recently proposed 
by us [lOj ] , determines the full set of critical exponents in 
terms of the fractal structure of open and closed physical 
lines. 
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